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A quantum chain model of many molecule motors is proposed as a mathematical physics theory 
on the microscopic modeling of classical force- velocity relation and tension transients of muscle fibre. 
We proposed quantum many-particle Hamiltonian to predict the force- velocity relation for the slow 
release of muscle fibre which has no empirical relation yet, it is much more complicate than hyperbolic 
relation. Using the same Hamiltonian, we predicted the mathematical force-velocity relation when 
the muscle is stimulated by alternative electric current. The discrepancy between input electric 
frequency and the muscle oscillation frequency has a physical understanding by Doppler effect in 
this quantum chain model. Further more, we apply quantum physics phenomena to explore the 
tension time course of cardiac muscle and insect flight muscle. Most of the experimental tension 
transients curves found their correspondence in the theoretical output of quantum two-level and 
three-level model. Mathematically modeling electric stimulus as photons exciting a quantum three- 
level particle reproduced most tension transient curves of water bug Lethocerus Maximus. 
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Muscle is a bundle of cylindrical muscle fibres. One 
muscle fibre enclose hundreds of thinner cylindrical fi- 
bre called fibril. A fibril is a train of sarcommeres con- 
nected head to tail. The sarcomere, as seen by light mi- 
croscopy, is a composite structure of actin filament and 
myosin filament(Fig. [I]). The two filaments are bridged 
by long myosin protein called myosin molecule motor (see 
Ref. or biology textbook for a more detail review). 
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FIG. 1: The longitudinal section of filament array in muscle 
fibre. Myosin molecules are long polymers as cross-bridges 
between the actin filament and myosin filaments. 

The length of a sarcomere is observed to decrease dur- 
ing muscle contraction. Huxley modeled the muscle con- 
traction as mutual movement between actin filament and 
myosin filament .4]. When the myosin molecule in the 
myosin filament is excited up to attach on the actin fil- 
ament, myosin molecule would combine with Adenosine 
Triphosphate molecules, and convert energy from Adeno- 
sine Triphosphate hydrolysis into mechanical force, this 
is the attached state. In the detached state, myosin 
molecule is at rest and doing nothing. The arrangement 
of myosin molecules along the length of filament is in- 
commensurate. A stochastic model was introduced to 
describe the cooperative behavior of molecules motors 
which has disordered arrangement along a backbone^. 
Experiment research is beyond the description of classi- 
cal mechanics two-state cross-bridge model 6] which as- 
sumed the probability of the two states — detached state 
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and attached state — satisfy a pair of coupled differential 
equations. 

Myosin molecules looks like a long arm ended by a 
head domain. The myosin in muscle cell is usually called 
myosin-II, for it has two heavy chains in the head domain. 
A second kind of myosin molecule with single heavy chain 
was found in non-muscle cells, they are termed as myosin- 
I. The average length of myosin-II is about 160 nm|2]. 
Visible light corresponds to a wavelength range of 400- 
700 nm. The electromagnetic wavelength comparable to 
the length of myosin-II falls in ultraviolet region. 

A muscle fibre surrounded by membrane is under con- 
trol of nerve cell. Membrane is a filter system which is 
permeable for certain ions but is impermeable for other 
ions. The imbalance distribution of ions across the mem- 
brane results in an electric potential difference of -60 to 
-90 mv. The electric signals from nerve cell can modify 
the permeability of membrane. An active muscle would 
generate electric signal. Scientist use this electric signal 
to represent the tension inside the muscle. If we insert a 
needle with two fine-wire electrodes into the muscle, the 
electric activity of muscle can be detected and recorded 
by an oscilloscope. This electromyogram is in widespread 
use for medical examination of muscle. 

Electromagnetic wave is termed as photon in quan- 
tum mechanics. The electric signal generated by ac- 
tive muscle is physically equivalent to a wave packet 
of photons. Quantum scattering between photons and 
molecules within biological system is a common phenom- 
ena, so does in muscle. Molecule only absorb photons at 
resonance frequency. We assume there exist a one-to-one 
correspondence between conformational change of mo- 
tor molecule and the hopping between quantum states of 
molecule. We define the excited quantum level as the at- 
tached states in which the motor molecule is extended to 
a longe shape. The detached state is assumed to be equiv- 
alent to a quantum states with lower energy. The motor 
molecule becomes shorter in detached states. Besides the 
two quantum states of molecules, another quantum state 
of photon is introduced to induce the hopping between 
the two quantum level. Photon can propagate among 
different molecules. It is photon that couples different 
motor molecules to work together. Quantum physics has 
many techniques to control quantum states, thus it is 
possible to control the conformation of molecule motors 
by electromagnetic wave. 

Muscle fibre can also be stimulated by chemical so- 
lutions which involves many complicate biochemistry 
There was an argument on chemical reactions in mus- 
cle using quantum mechanics Q, it is not related to my 
mathematical physics model. The time scale in my mod- 
eling is much larger than the chemical reaction. We only 
focus on the electrical stimulation of muscle fibre with- 
out adding any external chemical solutions. To make this 
mathematical physics model work for a real muscle fibre 
system, the following assumptions are necessary: (1) The 
strength of electric field should be so strong that it domi- 
nates the conformation of molecules over thermodynamic 



fluctuations. (2) For the muscle fibre without any exter- 
nal input of stimulus chemical ions, the conformational 
change of motor molecules governed by electromagnetic 
field is much larger than that induced by chemistry. This 
is because electric field is distributed across the whole 
space, it exerts a global force on the motor molecules. 
While chemical ion only interact with molecule locally, it 
does not come into action unless it finds the correct bind- 
ing site. (3) The spatial conformation of motor molecule 
can be modified by concentration of some chemical ions. 

We simplified the muscle fibre as a one dimensional 
chain of many giant quantum particles which represents 
molecule motors. The velocity of sliding actin filament 
acts as environmental parameter. These giant quantum 
particles are excited by absorbing photons, and relaxed 
by emitting photons. We arrange these particles reg- 
ularly along a one dimensional chain. The wavelength 
of the stimulus electromagnetic pulse is assumed to be 
much larger than a sarcomere. The configuration of spa- 
tial arrangement within a sarcomere is supposed to has 
no influence on the output physics. 

The article is organized as follows: 

In section II, we present the quantum Hamiltonian for 
deriving analogous force- velocity relation with respect to 
the coresponding macroscopic quantity. The empirical 
Hill's relation is consistent with the steady solution of 
Heisenberg equations. We predict the force velocity re- 
lation for the slow release case and the non-steady state 
solution. 

In section III, we study how the force decays after the 
stimulus of electric pulse. Exact solution of quantum two 
level model reproduced similar tension activation curves 
of cardiac muscle. The analytical solution of quantum 
three-level model can generate most similar curves of ten- 
sion transients of insect flight muscle. 

In section IV, the tension transients curve of skeletal 
muscle is mathematically reproduced from the point view 
of quantum coherent state. The experiment curve coin- 
cides with the third order projection of coherent state. 

In section V, we derived the self-coupled quantum 
Hamiltonian for solving strongly coupled differential 
equations. For this case, the sliding velocity depends 
on the internal states of the motor molecule. 

The last section is a brief summary and outlook. 



II. FORCE- VELOCITY RELATION DERIVED 
FROM QUANTUM HAMILTONIAN 

In muscle release experiment, the two ends of the mus- 
cle is fixed at its resting length. Electric pulses with 
different frequency is superimposed to develop tension in 
the muscle. In a quick release experiment, one end of the 
muscle is set free quickly. For a slow release, the length 
of muscle is shortening with one end of the muscle oscil- 
lating or moving under control with in a relatively longer 
time interval. The tension- velocity relation is the relation 
between two macroscopic quantity. In the following, we 
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shall express the force by microscopic density operator, 
while the velocity is taken as macroscopic parameter. We 
first assume the velocity of contraction is controlled by 
the density of some ions in the solution. We only model 
the excitation and decay process of muscle. 
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FIG. 2: A one dimensional chain of quantum two-level parti- 
cles. Each particle represents one myosin molecule motor. 



Force-velocity relation for quick release of 
muscle 



(d t + vd x )p l ab = ig a b(p. 



bb 



Paa)-^Pab ~ \P\b, (3) 



where Aw; = utf — oj^ is the energy difference between the 
two states. At time to, and position zo, the initial number 
of molecule motors in state | a) and | b) are denoted as 
rfi a and n b respectively, both rfi a and n b are functions of 
velocity v. In the steady states, the probability in both 
two sates are constants. Put the steady state constrain, 
dtPa a — and d t p l bb = 0, on the equations of motion — 
Eq. © and Eq. ([3]), one would arrive at an equality of 
probability distribution and velocity(detail calculation is 
presented in Appendix [Al . 



(52 PL 



N, 



)(v + sR(v)) 



N t s[n° - nl 



(4) 



We summarize a myosin molecule motor as a multi- 
level quantum particle. One eigenstate of the particle 
with respect to each level represents one conformation. 
In the excited state |o), myosin molecule is attached on 
actin filament. In the ground state the quantum par- 
ticle is at rest in myosin filament. One photon may excite 
the quantum particle from ground state to excited state. 
The particle hops from state \a) to state \b) by emitting 
one photon. The quantum representation of the mutual 
sliding movement between actin filament and myosin fil- 
ament is to annihilate one quantum state on one lattice 
site and generate it again at the next nearest neighboring 
lattice site. As the sliding motion has only one direction, 
switching two neighboring lattice sites would contributes 
a minus sign in front of the hopping operator. The quan- 
tum Hamiltonian of this one dimensional chain(Fig. [5]) 
is unitary, 

H = 52 U i\ a )i,i( a \ + u} i\ b )i,i( b \ +9ab\a)i,i(b\ + gba\b)i,i(a\ 



where R(v) is complicate function of velocity, 



a)i(a\ i+ i + \b)i(b\ i+ i - \a) i+1 (a\i - \b) i+ i(b\i].(l) 



where s is the unit lattice spacing, and v is the absolute 
sliding velocity. The hopping coefficient g a b{t) = p Q E(t) 
is electric field. Usually we choose symmetric coupling 
coefficients, g a b(t) — gba{t) = g- This Hamiltonian is nei- 
ther single-particle Hamiltonian nor tight-binding model. 
For the particle at lattice site i, the probability of being 
in the excited level |a), is measured by the density opera- 
tor, \a) it i{a\ = p l aa . p\ b = \b)ia(b\ defines the probability 



of being in state p l ah 



i t i{b\ is proportional to the 



complex dipole moment. The probability density func- 
tion are not independent, they must satisfy the constrain, 

Paa+Pbb = !■ 

The time evolution of the density operators are gov- 
erned by the equation of motion, p — — ih \H, p\. For 
the given Hamiltonian ([TJ, the equations are 

(dt + vd x )p\ a = n° a + igabplb-WbaPba ~ 2 ~P l aa, 

s 

(d t + Vd x )p\ h = n° b + igbaPba-WabPab ~ ^Pbb: ( 2 ) 
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where k — uj/c, c is the speed of light. /„ is the frequency 
of the stimulus electric pulses. N t is the total number of 
particles. We define the tension of the one dimensional 
chain as the total number of particles in excited states, 



Tension = x^^a' 



(6) 



where x is a renormalization coefficient, we set it as \ = 1 
for convenience. The myosin is attached on the actin fil- 
ament in excited states. The more myosin molecules are 
excited, the stronger physical bonds there would exist be- 
tween actin filament and myosin filament. The maximal 
tension is determined by the total number of particles, 
Tension max — N t , i.e., all particles are excited to con- 
nect the two filaments. 

The existence of a maximal tension value is consis- 
tent with experiment measurement. If one apply a single 
electric shock across membrane, the tension of the mus- 
cle fibre first increase from zero to a transient maximal 
point and then decays to zero. The electric shock is an 
energy packet with finite number of photons. Suppose 
the energy packet contains N p photons, and N p < N t . It 
can maximally excite N p molecule motors. In the mean- 
time, the motor is decaying and emitting photons. If 
Nt photons are sent into the muscle at one time before 
any decay process, all the molecules will be excited, then 
it will reach the maximal tension state. It always take 
certain time for the muscle to eat all the photons, more 
photons must be added to fill the loss of photon by decay. 

The relation between total probability density and ve- 
locity obtained for steady states, Eq. Q, is quite simi- 
lar to Hill's empirical force- velocity relation[7|. Here we 
know the physical origin of every parameter. Although 
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FIG. 3: The numerical plot of tension-velocity relation ac- 
cording to Eq. ((H) . 



there is a complex function R(v) involved in Eq. (j4|). 
the tension- velocity plots a hyperbolic curve, tension de- 
creases with increased velocity (Fig. [3]). The equations 
of motion, Eq. J2j) and Eq. (J3j) , suggest that the veloc- 
ity is actually a decay coefficient. Larger velocity leads to 
faster decay. Thus the number of physical bonds connect- 
ing the two filament decreases as velocity increases. This 
give us a natural explanation to experimental results. 



B. The equation of motion for the slow release 

Besides the decay induced by mutual sliding move- 
ment, there exist spontaneous decay due to thermal fluc- 
tuation and vacuum fluctuation. When we take the spon- 
taneous decay into account, the equation of motion for 
the density operator is given by 



P = -j\ H 'P\ 



1 



(7) 



where the T is friction tensor, I\j = jidij. The equation 
of motion has similar form as that for quick release, 



(S t + vd x )p\ a = n° a + igabPab-iQbaPia ~ j'aPla, 

(d t +vd x )pl b = nl + ig haP \ a ~ig ah p\ h ~i h p\ h , (8) 



(d t + vd x )p\ h = ig ab (pl b ~ pl a )~iAuj lP l ab - i abP l ab , (9) 

where Aw, = wf - u\, i ah = \{i a + j' b ). The decay 
coefficients of velocity now becomes the composition of 
velocity and spontaneous decay coefficients, 

7l = (2-+ 7 a), 7b = (2-+7fc)- (10) 
s s 

For the steady state, d t p l aa = , d t p\ b = 0, the relation 
equation between the excited state density operator and 
the composite decay coefficients is 
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where An(z, v, t) — n®j' a 1 — n®"f' b 1 is initial probability 
difference. R(v) is a complex function of velocity, 
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If the velocity term is much larger than spontaneous de- 
cay term, i.e., 2v/s ^ 7a, 7&, one neglects spontaneous 
decay, then Eq. (fTTj) reduced to former tension-velocity 
relation Eq. Therefore, Eq. @ is the force- velocity 
relation for the quick release of muscle. If the muscle 
is released slowly, the spontaneous decay may become 
equally important as velocity term, 2v/s 7a, 7b- In 
that case, the tension-velocity relation should obey Eq. 
(fTTj) . Thus Eq. (fTTj) is the general force- velocity relation 
for steady state. As for the non-steady state, we have 
to solve the equations of motion to find how the excited 
state density operator evolves as a function of velocity. 



C. The tension-velocity relation beyond steady 
state 

Insect flight muscle oscillate far more rapidly than the 
frequency of the input nervous impulse. A muscle is un- 
able to expand, it only contracts. The muscle expansion 
is achieved by contracting an opposing muscle. The os- 
cillation frequency of muscle is the frequency of the os- 
cillation between excited state and relaxed state. The 
wf|o)i,»(o| + uj\\b)i^(b\ terms in Hamiltonian Eq. (UJ) de- 
scribes two free harmonic oscillator which are oscillating 
with frequency u> a and w b respectively. The ultimate 
output frequency of Hamiltonian Eq. (UJ is regulated by 
the hopping term between the two levels and the mutual 
sliding term. We perform Fourier transformation on the 
density operator, 



W = 



iksl I I 

e (afc 



-iksl 



a fc >,(13) 



and substitute Eq. (fT3j) into Hamiltonian Eq. (jTj). The 
Hamiltonian now becomes more brief, 

H = H Q + J2(9ab\a k )(b k \+g ba \b k ){a k \), (14) 



where H is the exactly solvable part for regulated free 
oscillator, 
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(15) 



The regulated oscillation frequency depends on velocity 
and momentum vector k. The free oscillation modes has 
a maximal frequency and a minimum frequency, 
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The band width for existed free oscillation modes is de- 
termined by the sliding velocity. An oscillating electric 
field, 



9ab(t) = P E(t) = p E cos(/i), 



(17) 



would induced the hopping between two levels. The state 
vector for the Schrodinger wave equation, d t \^p(t)} — 
— ifir 1 H\ip(i)) , is written as 

2v 

\ipk(t)) = C a exp[-i(wfc sm(ks))t]\a k ) 

s 

+ C b exp[-i(u%- — sm(ks))t]\b k ), (18) 
s 

The density operator are given by P* a {t) — C*C a - 
Within rotating wave approximation, the terms that do 
not satisfy energy conservation are dropped, such as 
d + \a)i,i(b\ and d\b)i t i(a\. The eigenstate density operator 
under the first order perturbation of the hopping terms 
g a b{t) is 



p 2 E 2 sin"^ (Aa; - (2v/s) sin(fcs) - f)t] 
Ah 2 [Auj - (2v/s) sin(fcs) - f] 2 



(19) 



where Aoj = uj a — uj b . The non-steady state tension de- 
crease following v~ 2 regulated by an oscillation factor. 
The ultimate oscillation frequency of tension includes 
three sources, 



ujk = Auj — (2v/s) sin(fcs) — /. 



(20) 



When the frequency difference of the two level vanishes, 
Au> = 0, the output frequency is utterly dependent on 
the sum of perturbation frequency and sliding velocity. 
The fast oscillating ion flow of K + , Na + across the mem- 
brane cause the potential to rise to +40 mv and restore 
its original value within a few milliseconds @- This os- 
cillation frequency is coincide with some muscle oscilla- 
tion. Some insect flight muscle contract 1000 times a 
second. My conjecture is that the oscillation frequency 
of muscle has some possible relation with the oscillat- 
ing of ion flow. Certain insect flight muscle can oscillate 
more rapidly than the frequency of the input nervous 
impulse. Eq. ([20t told us the output frequency not only 
depends on the frequency of input nervous impulse, but 
also depends on the sliding velocity and the frequency 
difference between the two levels. The sliding velocity 
may increase the output oscillation frequency. The out- 
put tension measured in experiment is the superposition 
of different oscillator modes, 
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(21) 



The two ends of the one dimensional muscle fibre is al- 
ways attached to some tissues. If there are N+l particles 
along the chain, the beginning and the ending particle are 
fixed. As sin(fcs) in dispersion Eq. (|20[) is a periodical 
function, independent standing wave modes is confined 



in the domain (— it < ks < tt). The allowed wave modes 
are 

TT 2lT 37T TO7T (N t — 1)tT 
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(22) 



where L is the length of one sarcomere. All particles are 
oscillating in a standing wave pattern with frequency oj. 
For a selected modes the phase difference between 
nearest neighboring particles is If the sliding veloc- 

ity is zero, v = Q, the standing wave modes vanished. 



III. QUANTUM THEORY OF TENSION 
TRANSIENTS UNDER STRETCH ACTIVATION 

A stretch activation is to lengthen the muscle fibre by 
external force, or to release the fibre to certain length 
with in a extremely short time. The fibres were glued 
between two glass rods. One rod is attached to the anode 
pin of a force transducer. The length and tension signal 
are displayed on a double beam oscilloscope (Tl). 

Both the myosin filament and actin filament are long 
helical chains. The elementary structure to build myosin 
filament looks like golf club. They entwined around each 
other to produce a helix with diameter of about 15 nm. 
The actin filament is a double helix structure composed 
of globular protein. The local electric charge distribution 
along the myosin chain and actin chain are not homoge- 
neous. A mutual sliding between the myosin and actin 
chain would inevitable modify the local potential config- 
uration. The raising or lowering of local potential barrier 
induce the hopping of electrons among different energy 
level. Many photons are generated within very short 
time. The stretch activation by mutual sliding move- 
ment is equivalent to stimulating the muscle fibre with 
an electric shock pulse. 

The time evolution of tension induced by a stretch ac- 
tivation includes two periods. In the first period, there 
is a rapid decay due to the extremely large sliding ve- 
locity. A large number of photons are generated but are 
not absorbed by molecule motors. In the second period, 
the two filament stopped mutual sliding. The accumu- 
lated photon began coming into action. This period is 
dominated by the spontaneous emission and absorption 
of photon. 

The ideal definition for contraction velocity is 



v(tn) = lim 

At->0 



L(t + At) - L{t Q ) 
At 



(23) 



L is the length of sarcomere. In fact, what we measure 
in experiment is the average velocity within At. At the 
moment of quick release, the length experienced a sudden 
change, 



AL(t ) = / AL(t)6(t - t )dt. 



(24) 



The sliding velocity at t$ is almost infinity, v{fo) ~ oo. 
While after the stretch, the length does not change any- 
more, the sliding velocity becomes zero. The Heaviside 
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function provides us a good mathematical description to 
the length variation during a stretch. 

In the microscopic phenomena, a stretch activation 
would induce a sudden change of the particle distribu- 
tion on different energy levels. There exist an optimal 
distribution which acts as the minimal point of a poten- 
tial configuration. Any deviation from this optimal point 
has a tendency to go back. The optimal particle distri- 
bution is different for different types of muscle fibre. The 
external load applied upon the muscle is actually a high 
potential barrier. The particles are trapped in a local 
minimal point to counterbalance the external force. As 
a muscle is only able to contract, the intrinsic potential 
without external load is asymmetric. 



Quantum two level model for stretch activation 
of cardiac muscle 



The stretch activation of cardiac muscle has been ex- 
tensively studied in experiment 11]. Various different 
species of vertebrates were put into use, including Hum- 
ming bird, frog, Guinea pig, rat, rabbit. The tension 
transients following length perturbation were recorded 
by force transducer. At the instant lengthening the mus- 
cle, the force jumps up to a maximal force. Then the 
force drops to a temporary minimum. After several pe- 
riod of small oscillation, it finally reach a new maximal. 
The length of the muscle were kept invariant for a few 
seconds, and was suddenly release to its original length. 
A sharp drop of force appears immediately following by a 
fast recovery. The force transients for stretch and release 
looks more or less like mirror-image(Fig. 5. in Ref. [Tij). 
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FIG. 4: The mechanical oscillator as an analogy of the quan- 
tum two level model. 

A quantum two-level system provide a good model for 
the physics of force transients of cardiac muscle. A me- 
chanical analogy of this quantum model is an oscillator in 
thick liquid. As shown in Fig. [U a heavy ball is attached 
to one end of a elastic rod. The sudden lengthening of 
muscle is to hit the ball with a pulse. The ball pops up to 
hit the upper level. Since the length of muscle after the 
sudden stretch keep invariant, a spring from the upper 
level will attach the ball to prevent it from going back. 
However the elastic rod will pull the ball down, the com- 
petition between the elastic rod and the spring from the 
upper level induced a local oscillation. When the kinetic 
energy is finally consumed by the friction of the thick liq- 
uid, the ball stopped oscillating. For the sudden release 
of muscle, the pule hit ball in the opposite direction of 
the lengthening. As the molecule motor runs only in one 
direction, it is reasonable to introduce a bias potential. 
This bias potential makes the ball a little bit harder to 



restore from the bottom level. Here the bias potential is 
gravitational field. 

We assume the relaxed state of a muscle fibre is reached 
at an optimal ratio between attached states and detached 
states. One straightforward ratio is to set the number of 
particles in excited state and ground states as equal. This 
relaxed state is the vacuum state of the quantum two- 
level model. The vacuum state does not means there is 
no particle, but means the number of particle is exactly 
the number of antiparticle. In the mechanical analogy, 
this relaxed state corresponds to the equilibrium position 
of the ball without external pulses and gravitation field. 
If there exist bias field, the mechanical equilibrium is 
obtained by the balance between gravitational field and 
elastic force of the rod, nAh = mgAh. We assume the 
optimal ratio is determined by the chemical potential of 
the two states, 



P° 
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P° bb 



exp[ 
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(25) 



where /j, a and p b are the chemical potential of excited 
state and ground state respectively. In the following, 
the density matrix is actually the difference between the 
usual density operator and the vacuum density operator, 
we simply denote (p a /3 — p° a p) — > p a /3 for convenience. 
The vacuum density distribution determines the resting 
stiffness of a muscle fibre. 

A lengthening activation is equivalent to a positive 
pulse. All the particle in state \b) are driven to state 



AZf(to) = +Ni oa d, Paa(t ) = l, Pbb(t ) = 0. (26) 

The sudden release at t\ is negative pulse, all particles 
shift to state |6). 

AL{t 1 ) = -N loadl paa{h) = 0, p bb (h) = l. (27) 

Thus the length perturbation provide the initial condi- 
tion for time evolution of density operator. Since the 
length is constant after the stretch activation, we ignore 
the velocity term in Hamiltonian Eq. (JT]). The equation 
of motion including spontaneous decay reads, 

dtPaa = igabPab — igbaPba — JlPaa, 

dtPbb = igbaPba-igabPab+JlPaa, 

dtPab = igabiPbb - Paa)-iAu}p ab - J 2 Pab, 

dtPba = -igabiPbb - Paa)+iAojp ba - J 2 Pba, (28) 

where Au; = ui a — u . The coupling coefficient are chosen 
as the symmetric, g ab — g ba — g. We focus on the case 
of two degenerated states, i.e., Aw = 0. In the begin- 
ning, all particles are in state p afj (0) — & a b5p b - This 
corresponds to a quick release. Applying the standard 
ansatz 

4 

Pc/sit) = X] c «£ exp [ A 9*]> oc,P = a,b, (29) 

9=1 
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one reduces the differential equation to algebra 
equation [l2j]. It fi nally yields the time evolution of prob- 



ability in state \a) [l2j. 



Paa (t) 



2ff 2 



4.9 2 



7i72 



Aie Al * Are 



\ 2 i 



Ai — A2 A2 — Ai 



(30) 



This is the conventional solution in the textbook of quan- 
tum optics [12]. The plot of this solution for three dif- 
ferent cases were also shown in the book of Ref.[l2|]. I 
will reproduce some figures of Ref. [l2| in the following to 
show what the tension transients curve looks like. I will 
compare different cases with experimental curves which I 
have no copyright to use here. p aa {t) governs the tension 
transient. Here 



1 



1 



Ai,2 = -3(71+72) ± gVwi ~72) 2 - 16# 2 - 



(31) 



After the stretch activation, there is no more external 
stimulus. The molecules undergoes radiative damping 
through spontaneous emission. The time evolution of 
Paa(t) relies on the ratio between driven field g and damp- 
ing rate 7, here we have set 71 = 27, 72 = 7. If the 
driven field is much larger than damping field, there is an 
oscillation upon the exponential increase of tension(the 
upper curve in Fig. [5] (a)). In the opposite case, the 
damping dominates the dynamic evolution, it is an ex- 
ponential increasing without oscillation(the upper curve 
in Fig. [5] (b)). As the density operator must satisfy 

Pbb ~r Paa 

= 1, the time evolution of p bb = 1 - p aa just fol- 
low the image curve of p aa . For the intermediate case, the 
pumping strength is comparable with damping rate, both 
the lifetime and amplitude of oscillation are reduced, it 
approaches to an exponential increase for p aa (the upper 
curve in Fig. [5]) or an exponential decay for ^(the lower 
curve in Fig. [6]). 

Comparing these tension transients curve with the ex- 
periment records in Ref. [ll[, we find Paa(t) in Fig. 0(a)) 
fits well with the muscle tension transient under stretch 
activation, including the tension-time course of rabbit 
papillary muscle in Fig. 2 of Ref. [HI], the tension-time 
course of cardiac papillary muscle and insect flight mus- 
cle in Fig. 26 of Ref. [ilf . As for the release activation, 
the tension-time course of rabbit papillary muscle of Fig. 
26. A of Ref. [llj and Fig. 2 of Ref. [ll| fits very well 
with p aa {t) in the theoretical Fig. [5] 

The tension-time course of Humming bird muscle un- 
der stretch activation and release activation looks like 
mirror image, however the life time and amplitude of 
tension oscillation under stretch activation is larger than 
that for the release case(Fig. 5 of Ref. [H|). The stretch 
activation curve is well fitted by pbbit) — 1 in Fig. [5] 
(a)). The release activation curve fits better with p aa {t) 
in Fig. |51 The release activation is the tension recov- 
ery from state \b), the stretch activation is the recovery 
from I a). The asymmetric tension-time course indicates 
the damping rate for \b) is higher than the damping for 
state I a). In the mechanical analogical phenomena, the 



damping may comes from friction of liquid or the bias ex- 
ternal field. In polymer physics, the extension process of 
a long polymer is not always exactly the reverse process 
of contraction. I guess the motor molecule has similar 
behavior, this leads to the The asymmetric tension time 
course. We may call this a hysteresis phenomena. 




(a) 



(b) 



FIG. 5: (a) The upper curve is the tension transient of p aa for 
a strong pumping field and low damping rate, (7/7=3. The 
lower curve is the time evolution of {pbb — 1)- (b) The tension 
transition for a weak pumping field and high damping rate, 
g/7=0.3. The upper curve is p aa (t), the lower is {pbb — !)• 




FIG. 6: The time evolution of p aa {t) for 3/7=1. 5(the upper 
curve). The lower curve is (pbb(t) — 1) for 17/7=1. 5. 

The tensions time course for a weak pumping field and 
high damping rate(Fig. [5](b)), 5/7=0.3, reproduced the 
high tension state of rabbit psoas muscle. The tension 
transient of stretch and release are exactly mirror-image 
of exponential decay(Fig. 15 A of Ref. |ll|). 

The tension transient Eq. (f30|) has two extremal case, 
one is for extremely high friction 7 3> g, the tension first 
decreases, and then exponentially rise to a stable value. 
The other case is for very strong pumping field, g S> 
7. In this case, the tension shows swift oscillation, the 
envelop of the oscillation is an exponential decay curve. 
These two cases have no resemblance among the tension 
course of the experiment in Ref. [11]. There might exist 
some mechanism to prevent the cardiac muscle falling 
into these two extremal cases. 



B. Quantum three-level model for tension 
transients in insect flight muscle 

The tension transients of insect flight muscle under 
length perturbation is different from that of vertebrate 
including rabbit, humming birds, frog, pig, and so on. 
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FIG. 7: (a)The time evolution of p a a(t) for the case that 
friction 7 is much higher than pumping field g, 3/7=0. 05. (b) 
The time evolution of p aa (t) for g/7 1, the driving field is 
much higher than friction. 

The experiment in Ref. reported the muscle tension 
course of water bug Lethocerus Maximus. If the mus- 
cle is prepared in a regular solution without Ga ++ , the 
muscle tension almost has no response to length pertur- 
bation. In the presence of Ga ++ and Mg ++ , significant 
tension transient of Lethocerus Maximus muscle were 
observed in reaction to length perturbation [Tl| . 

The tension course curve under both stretch and re- 
lease activation showed a fast rise to a transient maximal 
and then followed by a slow decay [TTJ . The shape of the 
decaying tail is beyond the theoretical prediction of the 
quantum two level model. While a quantum three-level 
model fits the experiment data very well. This suggests 
that the molecule motors in insect flight muscle may be 
different from the motors in cardiac muscle. 

We still model the insect flight muscle as one dimen- 
sional chain of quantum molecule motors. Now the 
molecule motor has three states: the attached state |a), 
the metastable state \b) and the detached state |c). We 
assume the hopping between different states is induced by 
vacuum fluctuation. Since muscle works at room temper- 
ature, the stochastic environmental fluctuation may de- 
stroy the coherent correlation between different motors. 
The motors are almost independent from each other, thus 
the one particle model is a good approximation. The 
usual Hamiltonian for a three-level particle interacting 
with vacuum modes reads, 

H = ^9^ l(uia ^ k)t Ma}(c\+ g h 2 e^^ k)t d k \b}(c\] 

k 

-i $> f e - l( — fc)t 4 \c) H + g k 2 e~^~^d{ |c) , 
k 

(32) 



where uj a p = uj a — up is the frequency difference be- 
tween state \a) and \f3), a,f3 — a,b,c. d k and are 
the annihilation and creation operator of the kth vac- 
uum mode with frequency ui k . g k and g k denotes the 
coupling strength between vacuum modes and the three- 
level particle. 

The quantum three-level system has a similar mechan- 
ical analogy with the damping oscillator for the two-level 
model. The key difference is that the oscillator has a 
metastable partner suspended above the equilibrium po- 



— C- lb> 



1> 





FIG. 8: The mechanical oscillator analogy of the quantum 
three level model. 

sition(Fig. [5]). To perform a lengthening perturbation, 
we sent a shock pulse to hit the oscillator in the opposite 
direction of gravitation. The oscillator will rise to attach 
the partner ball and becomes a pair. The pair first rise to 
the attached state | a), and then decayed to the detached 
state |c). For a sudden release perturbation, the shock 
pulse hits the ball along the direction of gravitation field. 
The oscillator first drops to the detached state |c), and 
then is pulled up by the elastic rod. It continues to rise 
high enough to attach the metastable ball, and oscillates 
as a pair. This damping oscillator model gives us a rough 
phenomena of tension time course for insect flight muscle. 

We assume the myosin molecule motors as an equiv- 
alent three-level system. The tension of the one dimen- 
sional chain is proportional to the total number of par- 
ticles in the attached state, Tension oc p l aa . Another 
equivalent way of modeling this system is to take the 
whole chain as one giant three-level particle. The tension 
of the particle is proportional to the probability in the 
attached state, in this case, we take muscle as single par- 
ticle system which is relatively easy to handle with. An 
analytical investigation of Hamiltonian Eq. (|3"2"|) was well 
performed for the spontaneous emission of a three-level 
atom with quantum interference [13J. The dynamics of a 
quantum particle is governed by the Schrodinger equa- 
tion, dt\ip(t)) = —ih~ 1 H\ip(t)). The general state vector 
is written as 

|V(t))=C a |o)|0) + C 6 |6)|0)+5^Cfc4|c)|0). (33) 

k 

The molecule is initially in the attached state |o) and 
metastable state \b). When the frequency difference be- 
tween the attached state and metastable state is much 
smaller than the frequency difference between the at- 
tached state and detached state, i.e., ujba <C u>bc-,u ac , the 
Schrodinger equation admits an analytical solution 13], 



C a {t) 
C b (t) 



(^Aie Alt + BA 2 e A2t )e 



\2t\ p -(~f b /2+ioj ba )t 



{Ae Xlt + Be X2f 



-(7b/2)t 



(34) 



where the coefficients A and B are 

\ 2 C b (0) + 0.5^76<7 o (0) 



A = 
B = 



A2 — Ai 
AiC 6 (0) + 0.5V7a7&C7 a (0) 
Ai — A2 



(35) 



The decay rate of the two upper level are 7^ = 
2(7r ffl ) 2 L>(w fc ), 7a = 2(irg 2 ) 2 D(cj a ). D{uo) is the mode 
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density. Ai,2 determines the amplitude of the decaying 
metastable state, 



Al,: 



A7 



±\ ( 



A7 i0J ba , 



lalb 



(36) 



where A7 = 7& — 7 a is the difference of decay rates for the 
two upper level. The probability of being in the attached 
state is given by the density operator p aa (t). Here the 
weights of state vector, C a (t) and Cb(t), define p aa (t) and 
Pbb(t), 



Paa {t) = C* a (t)C a (t), p bb {t) = C* b (t)C b (t). 



(37) 



This density operator is actually a reduced density oper- 
ator by tracing out all vacuum modes, 



p aa = Tr k (\a, k)(a,k\), 



(38) 



To get a clear comparison with the two-level model, we 
map the Schrodinger equations of motion for C a [l3j into 
the equations of density operator, 



dtPaa = 

dtPbb = 

dtPab = 

dfPba = 



laPaa 



•\/7a76 / -tu),„, / 
1 v e 



P ba + e^Pab), 



Ibpbb 




la + 


7ft 


2 




7a + 


Oft 


2 



Pab 



Pba 



\Paa + Pbb), 



^ e \Paa + Pbb)- (39) 



The sum of the probability in the three levels satisfy the 
constrain, Paa +Pbb+Pcc = 1- For the two-level model, the 
sum of probability in the two level is conserved, dt[p a a + 
Pbb] = 0. While for the three-level model, p aa + Pbb = 
1 — pec, the sum of the probability in the upper two level 
depends on the probability in the detached state |c). The 
equations of motion (|39p have already incorporated the 
evolution of p cc . One may verify dt[p a a + Pbb] 7^ in Eq. 

The two excited states are coupled to the detached 
state by the same vacuum modes. The photon emitted 
by the transition from \b) to |c) may induce the hop- 
ping from |c) to I a). The photon may comes from one 
molecule, but absorbed by another molecule. Thus there 
exist strong quantum coherence between the two decay 
channel. A sudden lengthening or release of muscle leads 
to the mutual sliding between the charged filaments. A 
packet of photons are created immediately but are bot 
absorbed. The stretch activation and release activation 
set up the initial value of population at different level. 
The deviation from equilibrium position results in a tar- 
get potential barrier, which keeps a certain number of 
excited molecule from tunneling into lower level. The 
number of survived attached states determine the final 
tension value after decay. 

One qualitative conjecture of Ref. [ll[ is that Mg ++ 
is able to vary the Ga ++ sensitivity of the contractile 




(a) 



(b) 



FIG. 9: (a)Initial population in metastable state is p(,t,(0) = 
0.8. The attached state is /9 aa (0) = 0.2. Decay rate of 
metastable state is 75 = 7. Decay rate of attached state is 
7 a = 4. The frequency difference of the two excited states is 
Uba = 4. This curve corresponds to a stretch activation, (b) 
No population initially exist in the excited state p aa (0) = 0. 
Pbb(0) = 0.5, 7b = 7, 7a = 7, tJoa = 4. This corresponds to a 
release activation. 



a 

H 




(a) 



(b) 



FIG. 10: (a) Comparing with Fig. [9l the decay rates for 
attached state increases to 7a = 14. 75 = 7, u)b a = 4, pbb{0) = 
0.5, Paa(0) = 0.2. Tension decay becomes slower, (b) The 
decay rate 7a = 4. The frequency difference is reduced to 

OJba = 1. 76 = 7, Pbb{0) = 0.8, paa(0) = 0.2. 



system. When the Lethocerus Maximus muscle is in a 
regular solution with only Ga ++ ions, the tension reaches 
the delayed maximum within in 70 msecpd|. When 
Mg ++ is added, the time to reach the maximum ten- 
sion becomes much longer, it is about 500 msec[ll|- Fig. 
18 B in Ref. [ll| showed the tension time course in Ga ++ 
and Mg ++ poor contraction solution. The experiment is 
performed at temperature T=18°C with 5 mM ATP. The 
ion concentration is estimated as pGa=8~9. pMg=7^8, 
pH=6.5. Another Fig. 23 in Ref. ll| showed the case 
that the concentration of Ga ++ and Mq ++ is even much 

The tension 



11] 



higher than that for Fig. 18 B in Ref. 
time course of Fig. 23 in Ref. [11] is recorded at tem- 
perature T=18.5°C with 5 mM ATP, pH=6.5, pMg-5, 
pGa-6.8. 

My analogue theoretical conjecture is the damping rate 
of exited states is related to the relative ratio of Mg ++ 
to Ga ++ and their absolute value of concentration. In- 
creasing Mg ++ maybe reduce the decay rate. Higher 
concentration induces lower decay rate. This theoretical 
conjecture is based on the solution of three-level model 
in Eq. (|34p . We found two very similar tension time 
course to the experimental curve of Fig. 18 B in Ref. 
pl| . The first Fig. E](a) simulates the stretch activation, 
for which the initial probability in the attached state is 
Paa(0) = 0.2. The second Fig. [U(b) is for the release ac- 
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tivation, no particle exists initially in the attached state. 
To get the similar tension transient curve with respect to 
Fig. 23 in Ref. [ll[ from the three-level model, we raised 
the value of damping rate 7 a from 4 to 14 in the numer- 
ical calculation. The theoretical output of Fig. [10] can 
be identified as the experimental curve of Fig. 23 in Ref. 
11| . This theoretical curve in Fig. [TU] shows a slower 
decay process compared with the two decay curves lower 
concentration of Mg ++ and Ga ++ (Fig. [9j. 




(c) 



(d) 



FIG. 11: (a) The frequency difference between the metastable 
state and attached state cdba is very large, uiba = 14. pbb{0) = 
0.8, p aa (0) = 0.2, 7i, — 7, -y a = 4. (b) The frequency difference 
uiba increases to u ba = 40. p bb {0) = 0.8, p aa {0) = 0.2, 7 6 = 7, 
7a = 4. (c) LOba IS LO b a = 200. p b b(0) = 0.5, p aa (0) = 0.2, 76 = 
7, j a = 7. (d) uj ba is almost infinity, u> ba = 900. p b b(0) = 0.5, 

Paa(0) = 0.5, 7ft = 7, 7a = 7. 

The frequency difference uif,a indicates the oscillation 
between the two excited states. The exponential decay 
of tension is strongly regulated by this oscillation. For 
small value of uiba = 14, the tension has a dominant peak 
followed by small waves (Fig. [TT] (a)). If we ma gnify 
the tension course of experiment Fig. 19 in Ref. [Tl| . 
the theoretical Fig. [TT] (a) is coincide with experimental 
curve. This curve of Fig. 19 in Ref. 11] is obtained at 
temperature =23°C with 15 mM ATP, and pGa=8~9. 
pMg=7^8, pH=6.5. Comparing the experimental condi- 
tion with that for Fig. 18 in Ref. [ll|, one may see the 
temperature raise about 5°C, the concentration of ATP 
grows from 5 mM ATP to 15 mM ATP. My conjecture 
is that increasing ATP and temperature enhanced the 
oscillation between two excited states. 

If we make the w 6a larger, i.e., uiba = 40, the decay 
process shows more periods of oscillation(Fig. [TT] (b)). 
When ujba = 200, the tension-time course becomes a 
wavy curve of exponential decay(Fig. [UJ(c)). This wavy 
curve is quit similar to the tension transient of rabbit 
psoas muscle(Fig. 15 B in Ref. (HI). For an infinite large 
frequency difference, uiba = 900, the oscillation is com- 
pletely suppressed. The tension transients approaches 
to a smooth exponential decay(Fig. [UJ(d)), which looks 
very similar to the experimental curve in the high tension 



states(Fig. 23. C, in Ref. [TT|] ) . The high tension states 
of Lethocerus Maximus muscle is reached by increasing 
pGa from 8~9 to 6.8, and keeping pMg=7~8. From the 
theoretical point of view, a large frequency difference ujba 
produced a large energy barrier. The quantum interfer- 
ence between the two excited states vanished. In mind 
of the theoretical approximation we used, uiba <C uib c and 
ujba <C ui ac , a large uib a indicates that the excited states 
are almost decoupled from ground state. The dynamics 
of the three-level system is dominated by decay of excited 
state. Since no pumping field can pass an infinite barrier, 
the evolution of excited states is an exponential decay. 

When the frequency difference uiba vanishes, the two 
excited states becomes degenerated. Both the two states 
obey an exponential decay, 

la lb 

p bb (t) = A 2 + B 2 e~~ t{ - la+lb ^ + ABe~ t (~ la+lh ^ 2 . (40) 

The tension now is the sum of the probability of the two 
degenerated states, Tension — p aa (t) + Pbb(t), 



Tension — A 2 ( — 

la 



l) + B 2 (^ 

lb 



l)e 



-*(7a+7b) 



(41) 



The tension evolution is an exponential decay plus a con- 
stant term. The constant tension depends on the initial 
value of the density operator. For this special case, three- 
level model with two degenerated excited states gives 
similar result as the two-level model. However, there 
is always an induced decay rate by emitting photons in 
three-level model due to the metastable state. In two- 
level model, the emitting photons can be brought back 
to excited state without violating selection rules, thus 
the output of two-level model does not produce all of the 
decay curves as the three-level model did. 



IV. PROJECTING QUANTUM COHERENT 
STATE FOR SKELETAL MUSCLE FIBRE 
UNDER ELECTRIC STIMULUS 

Comparing with cardiac muscle and insect flight mus- 
cle, skeletal muscle has two special characters: one is 
low resting stiffness, the other is its ability of extending 
much beyond its optimal length. When we summarize 
these two special characters into theoretical modeling, 
low stiffness indicates there are less number of molecules 
motors in the attached state when the muscle is not ex- 
cited. Two filaments connected by less bonding is easy to 
extend. The electrical signal representing muscle tension 
is physically explained as polarized photon. The strength 
of electric signal is proportional to the density of photon. 
Quantum coherent photons exist in many living biolog- 
ical system [ljj]. Even though an isolated muscle fibre is 
not alive, we can still shed some new light on the twitch 
tension evolution in reaction to the stimulus of electric 
signal from the point of view of coherent states. 
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External stimulus field will generate tension in the 
muscle. If the two ends of a muscle fibre is fixed, the 
tension developed by external stimulus is called isomet- 
ric tension. The isometric tension of skeletal muscle was 
measured many decades ago@. A superamaximal rect- 
angular pulses is used to stimulate the muscle. Both 
the second and third derivatives of twitch tension are 
recorded by some software @. 

We take the same quantum Hamiltonian Eq. (fT4)) of 
one dimensional chain of quantum two-level particles for 
modeling skeletal muscle, 

H = H + Y,(9ab\a k )(b k \ +9ba\b k )(a k \)- (42) 



The exactly solvable Hamiltonian Hq is given by Eq. 
(fT5"j) . A quantized electric field motivated the hopping 
between the two levels, 



g ab (t)=p E(t)=p [d(t) + d+(t)}, 



(43) 



where d is annihilation operator of photons. Applying 
the conventional rotating wave approximation, the total 
Hamiltonian reads, 



H = H +p Q Pabd + p o p* ab d + , 



(44) 



where 



Pab 



y 

k 



/ ; 1 a.fc){6fc j , Pba 



k 



\b k )(a k \. (45) 



Hamiltonian Eq. (|44p shares the very similar formulation 
with the general Hamiltonian for a coherent state[l5j]. If 
there is no special projection to certain number of pho- 
ton, we would always have a coherent state for photons. 
However, one must project the coherent state to a fixed 
number of photons to get similar curve to the tension- 
time of skeletal muscle. We denote the off-diagonal den- 
sity operator p ab by a brief symbol p — \pab- The 
Schrodinger equation, 

\ip(t + At)) = \ip(t)) - ih^HAt^t)), (46) 

may be expressed by Fock states, 

.At, 



\n,t + At) = \n,t) - i—[p\/n\n-l,t) - p*Vn + l\n,t)]. 

Here we have input the annihilation and creation rule of 
photon operator, 

d+ | n) = Vn+l | n+1), d | n) = s/n | n - 1). (48) 

Eq. (|4"7) is equivalent to a radioactive decay equa- 
tion from which one can get the standard Poissonian 
distribution [8]. Physicist define coherent state as the 
eigenstate of annihilation operator d. The coherent state 
with respect to Hamiltonian Eq. (|4"4")) is \pt), 

\pt) = exp[p*td + - ptd]\0), d\pt) = pt\pt). (49) 



The general expression of wave function as solution of 
Schrodinger equation is \ip{t)) — exp[^H t]\pt) . Project- 
ing the coherent state \pt) to Fock space gives us Posso- 
nian distribution, 



\pt) 



1 



{ P tr 



cx V[ --p*pt 2 ]j2 . , 

2 n=0 



n). 



(50) 



This coherent state is the superposition of many pho- 
ton states. A photon is the electromagnetic radiation 
of quantum particles. Usually the electromagnetic ra- 
diation of human muscle is in the infrared region at a 
wavelength of tens of microns. The detector of electric 
signal in the muscle is metallic needle. We assume the 
detector of photon has an optimal window which only 
allow finite number of photon pass, it can not capture 
all the photons. Then the projection of coherent state 
to this optimal photon number state is what the exper- 
iment measured. This electric signal is proportional to 
the amplitude of the optimal photon state, 

T(n,t) = \(n\pt)\ = {P > exp[-H 2 /2]. (51) 

V(« - n j! 

The experiment of Ref. Q only reported one complete 
tension-time course(Fig. 1 E in Ref. @). The tension 
transient of Fig. 1 E in Ref. Q can be theoretically fitted 
by 



T(6,r) 



(par) 3 



exp[— |par|], ar = t 2 



(52) 



where a is an adjustable parameter to rescale the unit 
time scale. One can choose different a to get different 
amplitude for the curve. Both the second and the third 
derivative of T(6, ar) looks very similar to experiment 
curve(Fig. [T2l . except the minor oscillating tail of the 
curve. This oscillating tail is most likely due to fluctu- 
ating experimental environment or the poor quality of 
equipment in 1968. Because the experimental curves in 
1977 of Ref. [ll| rarely showed similar wavy tails. 




FIG. 12: The tension transients T. The first derivative of 
tension is V — d T T. The second derivative of tension is A = 
d 2 T T. All the three curves are quit similar to those in Fig. 1 
E of Ref. 0. 



12 



V. THE SELF-COUPLED HAMILTONIAN FOR 
THE MICROSCOPIC REPRESENTATION OF 
VELOCITY 

There exist relative sliding motion between two fila- 
ments. The total number of particles described by the 
quantum Hamiltonian is sandwiched between the two fil- 
ament. For an almost infinitely large system, the sliding 
motion of filament acts as environmental parameters. So 
we summarize the sliding velocity as a constant parame- 
ter in the model. In the quantum Hamiltonian of the one 
dimensional chain model, the total probability of being in 
the two quantum states is conserved. The quantum den- 
sity operator measures the probability for a particle being 
in certain state. For the two-level model, the density op- 
erator follow p aa + pm = 1. For the three-level model, the 
density operators of three levels obey p aa + p^ + p cc = 1 . 

When the one dimensional quantum chain has finite 
length, the sliding velocity of filament may depend on 
the internal states of motor molecule that drags the fil- 
ament. The total number of particles sandwiched be- 
tween the two filaments is proportional to the length of 
the overlapping region. When the myosin molecules be- 
tween the two filaments are excited, they becomes active 
to drag the filament. However if the molecules falls out 
of the overlap region, even if they are excited, there is 
no active bond connecting the two filaments. We call the 
motor molecules out of the overlap region as unemployed 
molecules. The molecules within the overlap region are 
named after employed molecules. The mutual sliding 
movement make the quantum chain an open system. The 
total number of employed particles are increasing during 
contraction. 

A faster sliding velocity means the total number of 
employed particle increase or decrease faster, so we de- 
fine the mutual sliding velocity as the speed of increasing 
employed molecules plus a initial velocity v c , 

dN e . . 

v(t) = -^ + v c . (53) 

The absolute value of v(t) is equivalent to the decrease 
of unemployed molecules. We assume the initial velocity 
v c is constant. The initial velocity maybe induced by 
chemical or biological effect. For the simplest case, one 
may take v c = 0. The time dependent sliding velocity for 
the quantum chain of two- level particles can be expressed 
in terms of density operator, 

v = d t ^2(p 3 aa + p{ b ) + V c , p J aa + p{ h ± const. (54) 

3 

According to Heisenberg equation, p = — i[H, p], this ve- 
locity is given by the commutator equation of Hamilto- 
nian operator, 

v = -iJ2lH, f>U 4b] + «c (55) 

3 3 

Usually a quantum Hamiltonian is expressed by a func- 
tional of many operators which does not include the 



Hamiltonian operator itself. Here the expanding se- 
quence of operators for Hamiltonian Eq. ([1]) includes 
a velocity which is a functional of Hamiltonian operator 
itself. The explicit formulation of Hamiltonian Eq. ([1]) 
including the Heisenberg equation of velocity reads, 

H = ^2^\a) hl (a\ +(4\b)i,i(b\ + 9ab\a}i,i(b\ + 9ba\b)i,i(a\ 

+ -(1 - \a)i{a\i)H(\a)i{a\ i+ i - \a) i+ i{a\i) 
s 

+ 1(1- |6) i <6| i )F(|6) i <6|m -\b) i+1 (b\i), 

+ -^[l a )i( a li+i + \b)i(b\i+i - \a) l+ i(a\i - \b) i+1 {b\i]. 

We may call this Hamiltonian a self-coupled Hamilto- 
nian. Transforming the self-coupled Hamiltonian into a 
conventional Hamiltonian must follow the algebra rules 
of operating matrix. We denote the projection of Hamil- 
tonian to local eigen states as, 

(a\iH\a)i = wj 1 , (a\iH\a) l+1 = f? i+1 , 
{b\iH\b)i=wl (b\iH\b) i+1 =t b iii+1 . (56) 

tf i+1 and t b i i+1 is the hopping rate between two neigh- 
boring sites, uif and <jj\ is the local eigen energy of the 
particle at the ith site. After some algebra, we obtained 
the conventional expansion of the self-coupled Hamilto- 
nian, 

i 

+ [^ + tl i+1 s- 1 ]\b) i , i {b\ 

+ 9ab\a)i,i(b\ + g ba \b} l ,i(a\ + [— l -]\b) l (b\ i+1 

s s 

- ^[|a) m (4 t + \b)i + i{b\i] + - ^HlahH+r, 
s s s 

where G is the renormalization operator, 

q = s- X)[l a ^ a l»+i + \b)Mi+i ~ \a) i+1 {a\i - |6><+i<&|<]. 

i 

We can map the self-decoupled Hamiltonian expressed by 
density operator into a more familiar equation for quan- 
tum field theory Each of the two quantum states is rep- 
resented by a boson operator, 

(on\ = Q!j, \a)i=af, a = a,b. (57) 

Notice here we study an open finite system, the total par- 
ticle number is not conserved, i.e., afai + bfbi ^ const. 
The decoupled Hamiltonian is simplified into a two- 
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component boson system, 



,iv c ui: f. , ,iv c a;,-., 

+ \-f - -f K a * +i + i-f - -f\ h i^ 

+ gabaibi + g ba bta t - -y [a^a, + bj +1 b t ], 

1 - h E>+«m + b t b i+i - a t+i^ - b i+M ^ 

This nonlinear Hamiltonian has no resemblance in con- 
ventional quantum field theory due to the renormaliza- 
tion operator G. If we assume 1/s <C 1, then the operator 
G is approximated by a constant G ~ 1. Such a approx- 
imation is not always appropriate, as the hopping step 
size s is usually very small. For the other extreme case 
of 1/s ^> 1, the constant term 1 in the nominator of G 
can be neglected, 



G 



afai+i + bjb l+ i - a+ +1 a t - bf +1 bi ' 



(59) 



The key rules for handling this inverse quantum opera- 
tor G is the elementary commutator relation between a 
matrix operator and the inverse of another operator, 



-H — Cij. 



(60) 



We can assume if i ^ j, C^ — 0, this is a physical 
assumption since two elementary bosonic operators are 
sitting at different lattice sites. While for i = j, Cij de- 
pends on the specific formulation of the mattix d, and 
6j. We know only one universal relation, fij^- = 1, i.e., 
J-] = 0. Physically speaking, fij annihilate one par- 
ticle at the ith site. Its dual operator -i- is not the con- 
jugate operator a\. More over, it is unknown how to 
decompose the commutator between one elementary op- 
erator and the inverse of many composite operators, such 
as, 



-} = C(a,b). 



(61) 



Maybe we should first find out how to map the inverse 
operator into a series expansion of the conventional poly- 
nomial operators. 

In fact, quantum operator provides a natural descrip- 
tion of directed molecule motors. We can define a walking 
operator w^ 1 , when it operates on a spatial coordinate, it 
walks one step forward or backward, 



w + \x) = \x + s ), 



w 



\x) 



so)- (62) 



s is the unit step size. A directed molecule motor is 
characterized by only annihilation operator or creation 
operator, 

(w+) m \x) = |x + ms ), (w~) m \x) = \x-ms ). (63) 



The velocity of the motor is denoted as the number of 
step m within one unit time interval. For the sliding 
filaments of muscle fibre, m has a upper bound and lower 
bound. This is an open system in which the Hamiltonian 
is only consists of either annihilation operator or creation 
operator. This open physical system is strongly coupled 
with another chemical system or biological system. 



VI. CONCLUSION 

When single muscle fibre is teased away from the whole 
muscle, it is no longer a living organism. In the eye of a 
physicist, one muscle fibre is a long cylindrical material 
with complex internal structure. The polymers, chiral 
proteins and ions inside the fibre are properly organized 
both in space and functionality. An input electric signal 
induced the deformation of molecule motors, which in 
turn shortens the fibre. Generally I believe the electro- 
magnetic interaction(or in other words chemical bond- 
ing) determines the equilibrium conformation of myosin 
molecules. The stochastic fluctuation from thermal en- 
vironment only results in some minor modifications of 
local structure, even though sometimes this minor modi- 
fication plays a crucial role. One myosin motor molecule 
can be simplified as a giant quantum particle with finite 
quantum states. Each state corresponds to one stable 
molecule configuration. The input electric field induce 
the hopping among different states. 

We mathematically model a muscle fibre as a one 
dimensional chain of quantum particles representing 
myosin motor molecules. Besides the two quantum 
states, exited state and ground state in analogy of at- 
tached state and detached state of classical model, an 
external quantum state of photon is introduced to cou- 
ple the dynamics of one quantum particle with another. 
The mutual sliding between two filaments is included in 
the directed hopping operator of a quantum Hamilto- 
nian. Quantum dynamics of this quantum Hamiltonian 
system gives us similar force- velocity relation to the clas- 
sical model for the quick release case. The force- velocity 
relation for the case of slow release and unstable states 
is beyond the capable description of classical mechanics 
model. However using this quantum Hamiltonian, we 
can derived a rather brief force-velocity relation for the 
case of slow release and unstable states. The quantum 
equation of motion suggests that the mutual sliding is a 
source of the decay of excited states. Another source of 
decay is vacuum fluctuation or thermal fluctuation, we 
summarize them into a spontaneous decay process. The 
mutual sliding movement means breaking bonds and es- 
tablishing new bonds between two filaments. The ra- 
tio between the life time of the sliding induced decay 
and spontaneous decay determines the behavior of force- 
velocity relation. For a quick release, the sliding decay is 
much larger than spontaneous decay. For a slow release, 
the spontaneous decay dominates. Both the experimen- 
tal measurement of quick release and slow release only 
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records the behavior of steady states solution. For the 
non-steady state, the force-velocity relation is an oscil- 
lating function of time. There is a difference between the 
frequency of force oscillation and the input frequency of 
electric stimulus. This theoretical conclusion is consis- 
tent with biological observation — some insect flight mus- 
cle indeed can oscillate more rapidly than the frequency 
of input nervous impulse. 

We proposed the mathematical definition of the ten- 
sion of a muscle fibre as the total number of particles in 
the excited states. Thus tension transients is described 
as the evolution of the density operator of the excited 
states. Based on this hypothesis, we apply the conven- 
tional quantum two-level model and three-level model to 
check the tension transients course of muscle fibre. The 
quantum two-level model reproduced most of tension- 
time courses of cardiac muscle under stretch activation. 
The three-level model can theoretically generate most 
curves of the tension transient of insect flight muscle. The 
quantum three-level model also produced some tension- 
time course of rabbit psoas muscle for which the two- 
level model can not account. The key difference between 
three-level model and two-level model is that three-level 
model takes into account of the induced decay due to 
metastable state while the two-level model does not. The 
quantum three-level model can also generate the tension 
transients of a skeletal muscle fibre under stimulus of elec- 
tric pulse. However we checked the tension transients of 
skeletal muscle from a different point of view of quan- 
tum coherent state. The tension, the first derivative of 
tension and the second derivative of tension are all co- 
incide with the third order projection of coherent states, 
ft is not fully understood in physics why it projects to 
the third order, maybe there are better ways to model 
skeletal muscle. 

In the muscle stimulus experiment almost 35 years 
ago[llj], the tension-time course varies according to the 
concentration of Ga ++ and Mg ++ . In the quantum chain 
model, those mathematical parameters can varies to pro- 
duce similar curves to the experiment data. If there are 
too many theoretical parameters, it is big challenge to 
understand the physical or chemical relation between the 
concentration of Ga ++ and Mg ++ and those parame- 
ters. Fortunately there are only three effective param- 
eters in the quantum model related to the muscle ten- 
sion transient : the two decay rates of excited states and 
the frequency difference between the two excited states. 
The decay rates controls the long time behavior of ten- 
sion transient. The frequency difference introduces wavy 
oscillations upon the conventional exponential decay or 
exponential increase. It needs a series of systematic ex- 
periments to find out how ions control the three param- 
eters, the existed experiment data is not enough. One 
experiment fact is that if solution does not contain Ga ++ 
ions, delayed tension changed disappeared, if the Ga ++ 
and Mg ++ are added, distinctive tension transients were 
observed [llj- It seems as if Ga ++ ion drive the myosin 
motor molecule into metastable conformation and trans- 



form the molecule into a three-level quantum particle. 
While Mg ++ seems like to fix the conformation of myosin 
molecules so that they are invariant under stretch acti- 
vation. Maybe the ratio between the concentration of 
Ga ++ and Mg ++ is a possible controller of the two de- 
cay rate of the two excited states. Both theoretical and 
experimental works are required to solve the puzzle. 
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VIII. APPENDIX 

Appendix A: The derivation of force velocity 
relation for the case of quick release 

I present the detail calculations of force- velocity equa- 
tion here. The assumptions and approximations used 
in calculation are explained at every step. I first show 
how I derived the differential equations from the quan- 
tum chain model. Then I add the general procedure for 
solving similar differential equations in Ref. [l(| for the 
benefit of readership. The assumption used in solving 
the very similar differential equations have more detail 
explanation in laser physics theory [Toj. 

For a quantum particle, the probability of being in 
a quantum state is given by the square of the norm of 
\a), i.e., p —\ | a) | 2 . This probability is well defined by 
density operator, p aa = \a)(a\. For a two level system, 
there are four density operators, 

PL = HmH Pbb = \b)i,i{ b l 
Pab = \a)u(b\, p\ a = \b)i,i{a\. 

(Al) 

The diagonal density operator satisfy the constrain, 

l = i[|o) M <o| + |ft) M <6|], (A2) 

The equation of motion for these density operator is 

P=-j[H, P ] = - l -[Hp- P H]. (A3) 

The Hamiltonian of the one dimensional chain model is 
H = Y, 0J t\a)i, l (a\+uJ b i \b) h i(b\+g ab \a)^(b\+g ba \b) itl (a\ 

1 

11) 

+ —[\a)i(a\ i+ i + |6)i(6|i+i - \a) l+1 (a\i - |6) i+1 (6|j]. 
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Substituting this Hamiltonian and density operator into 
Eq. (|A3|i . one may get the equation of motion of the 
probability density in different states. Here we take p l aa 
as an example, 

tv 

ihPaa = — - \a)i](a\i + \a)i[(a\ i+ i - (a\i]} 

ID 

+ 9ba\b)i, l {a\~ g a b\a) l .i{b\-2 — \a)i{a\ i . (A4) 

s 

If the lattice spacing s is very small, we summarize the 
first term of Eq. (|A4|) into a continuum limit, 

-{[|a)i+i - |a)i](a|» + |o)i[(o|i + i - (a\i]} = d x p aa (A5) 



Here h is set to h = 1 for convenience. Then Eq. (IA4 
becomes more compact, 



(d t + vd x )p l aa = n° a + igabPU-WbaPba - 2-pl a (A6) 



Following a similar calculation, we get the equation of 
Pbb and Pab, 

(d t + vd x )p l bb = n° b + igbaPla-WabPlb - 2-P6 6 ,(A7) 



(d t + Vd x )p\ b = ig a b(Pbb ~ P l aa)~ i ^iP l ab dra ! t - 2 -p l ab, 

S (A8) 

where Aw; = (w^wf). Since p* ab = p ba . The equation of 
motion for p ba is just the conjugate of p ab . Mathematical 
theory told us, if the equation bear the form of 



, d d , . 

dt +v Jfa)f( x > t > v ) = 9(x,t,v), 



(A9) 



the solution is formulated as 



f(x,t,v) 



dt'g(x',t',v), x' = x-v(t-t'). (A10) 



The above equations of motion for density operator has 
exactly this kind of form. The solutions can be obtained 
by analogy of the calculation in quantum optics [10|. 
Within the rotating wave approximation, the coupling 
coefficient are 

9ab(t') = y E En^e-^'+^Ujz'), (All) 



Integration of Eq. 



gives out 



pi b (z,v,t)=i / dt'e*p[-(iu + 2-)(t-t')] 



xg ab (z\ t')(pl b (z', v, if) - pl a (z', v, t')). (A12) 

Following the assumptions of laser theory [Icj, the pop- 
ulation difference p bb (z', v, t') — p i aa (z' ,V,t'), the electric 



field mode amplitude E n and the phase <j) n vary little in 
the time s/(2v), then p ab is rewritten as, 



Pab = ~\^ E n exp[-*(/n*' + 4>n)} 



2«, 



dt'U n (z') eacp[-[*(w - f n ) + — ](t - t') 



x[pl b (z',v,t')-pi a (z',v,t')} 



(A13) 



U n (z') is the normal- mode function in space. For the 
one dimensional chain model, U n (z') is the distribution 
of electric field along the length direction of the chain. 
As the electric field is oscillating in space, one usually 
take sinusoidal z dependence: U n (z') = sin[fc n z']. As all 
know, exp[ik n z] — cos[k n z]±i sin[k n z], U n (z') is decom- 
posed into the sum of two opposite waves. For a sym- 
metric velocity distribution W(—v) = W(v), the quick 
cos[/c„z] terms are canceled in the integration, it finally 
yields U%(z) — 1/2. The density operator p ab now reads, 

p\ b = -iP±E n e- i ^ t +^(pi b -pi a )(D + +D_), (A14) 

where D± is the propagator function, 

1 



D± = 



2v/s + i(uj — f n zLkv) ' 



(A15) 



where k = lu/c, c is the speed of light. Substituting Eq. 
(IA14P into the equation of motion of p aa and p bb , one 
obtain the rate equation, 

dtp a a = K- R(v){p l aa - pl b )-2-p l aa , 

s 

dtp\ b - nl + R(v)(p i aa -p\ b )-2-p\ b . (A16) 

s 



where 



R(v) 



[G- + G+] 



8 h? 2v [ 

G± has the similar form of Lorentzian function 

4v 2 /s 2 



G ± 



Av 2 /s 2 + (u- f n ±kv) 2 



(A17) 



(A18) 



For the steady state, d t p aa — 0, d t p l bb — 0, the pair of 
rate equation leads to 



Paa Pbb 



An(z) 

TTrI 



(A19) 



Here An(z, v, t) is the initial population difference in the 
absence of field oscillation, 

An(z, v,t) = £- [n° a (z, v, t) - n° b {z, v, t)} . (A20) 



Eq. (|A19|) gives the relation between population differ- 
ence and decay coefficient. In our one dimensional chain 
model, the decay coefficient is just the linear function 
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of velocity. So Eq. (IA19[) is actually the force-velocity 
relation, 

{Paa- Pbb)(v + Rs) = . (A21) 



The rate function R(v) is a complex function of veloc- 
ity. At the resonance points lu — f n = ±kv, R(v) gained 
maximal value. 
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